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We study the electronic and magnetic structures of multiferroic CaMn 70 i 2 by first-principle 
calculations, based on the experimentally determined modulated crystal structure. We confirm 
the presence of a 3d orbital modulation of the Mn2 (Mn^^) sites, previously inferred from the 
Jahn-Teller crystal distortions. Our results indicate that in the multiferroic phase the magnetic 
structure of the Mn3 (Mn^^) sites is anharmonically modulated via orbitally-mediated coupling with 
the structural modulation, and that the Dzyaloshinskii-Moriya and exchange striction mechanisms 
contribute equally to the polarization. 

PACS numbers: 75.85.+t, 77.80.-e, 


I. INTRODUCTION 

Research into magnetoferroelectrics, in which im¬ 
proper ferroelectricity is induced by magnetic ordering, 
has had a significant impact in condensed mattei^hysics 
and materials science over the last ten years .^1® These 
materials have attracted wide attention because of their 
novel physics and their potential applications in multi¬ 
functional devices. One important and currently limit¬ 
ing factor is that, in order for these materials to be use¬ 
ful in technology, the magnitude of the electric polariza¬ 
tion must be of the order of that observed in proper fer- 
roelectrics. The multiferroic GaMn 70 i 2 (labelled GMO 
in the remainder) was shown to exhibit giant magneto- 
ferroelectricity of magnitude ^2.9 mG/m^ [IZHS] — one 
of the largest measured values of magnetically induced 
electric polarization. 

GMO has a trigonal strucuture (space group R3) be¬ 
low ^ 400 K, consisting of three symmetry inequiva¬ 
lent manganese sites that we label Mnl, Mn2, and Mn3 
— see Fig. [^a).l^ On decreasing the temperature be¬ 
low Too=250 K, the average crystal structure adopts 
an incommensurate modulation (IGM structure) prop¬ 
agating along the hexagonal c axis with structural prop¬ 
agation vector = (0,0,2.077) at 150 K. On fur¬ 
ther cooling, long range magnetic order develops be¬ 
low Tni ^90 K with a magnetic propagation vector 
Qm! = (0, 0,1.037)(AFM1 phase), followed by a second 
magnetic phase (AFM2 phase) below Tn 2 , described by 
(at least) two propagation vectors q ^2 = (0, 0, 0.96) and 
QmS = (0,0,1.12). Neutron powder diffraction showed 
that in the AFMl phase, all spins he in the ab plane 
with the two Mn^+ sites (Mnl on the A' perovskite site 
and Mn2 on the B site) forming ferromagnetic triangular 
layers that are rotated by 124° with respect to neighbour¬ 
ing layers along the c axis, as described by q^i- Spins 
on the Mn^+ (Mn3) sublattice are roughly antiparallel to 
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those in neighbouring Mn2 sites, but the precise angle of 
rotation is poorly constrained by experiments. 

Recently, it was proposed that in AFMl, the structural 
and magnetic modulations lock in a magneto-orbital he¬ 
lix and that the incommensurate orbital modulation, in¬ 
ferred from the Jahn-Teller distortions, is crucial in sta¬ 
bilizing the chiral magnetic structure.!^ In this model, the 
giant magnetically induced electric polarization would 
arise primarily from the inverse Dzyaloshinskii-Moriya 
(DM) effect .im This picture seems to be contradicted by 
first-principles calculations, which were performed on an 
average, unmodulated structure of GMO, and led to the 
proposal that the spin directions of Mn3 are stabilized 
by strong DM interactional^ rather than orbital modula¬ 
tions. Given this spin structure of Mn3, the electric po- 
larizatio n wou ld then be induced by symmetric exchange 
striction. b^F^ ljj^ paper, we revisit the magnetic struc¬ 
ture and magnetoelectric coupling mechanism in GMO 
using first-principles calculations on a series of structure 
models which approximate locally the modulated struc¬ 
ture. Our calculations confirm the existence of an or¬ 
bital modulation as described in our previous paper.l^ 
Furthermore, we propose that the magnetic structure of 
Mn3 should be modulated anharmonically with two wave 
vectors q^i and 3qmi- Finally, we show that the giant 
magneto-ferroelectricity observed in experiments is due 
to equal contributions from both exchange striction and 
the inverse DM mechanism. 


II. METHODS 

Our first-principles calculations are based on the 
density-functional theory impleme nted in the Vienna ab 
initio simulations package (VASP) We use the spin- 
polarized generalized gradient approximation with on¬ 
site Goulomb interactions U included for Mn 3(i orbitals 
(GGA+C/)P The projector augmented-wave (PAw;P 
method with a 500 eV plane-wave cutoff is employed. A 
2x2x4 Brillouin zone mesh is used. Electric polarization 
is calculated using the Berry phase method.^ 





2 


The magnetic structures of the Mnl and Mn2 sub¬ 
lattices in the AFMl phase, which have been well deter¬ 
mined by experiment, form the basis for our calculations. 
The magnetic propagation vector of the AFMl phase 
is approximated by the commensurate vector = 

(0,0,1). This approximation results in a 120° rotation 
(rather than the experimental 124°) between neighbour¬ 
ing Mnl and Mn2 layers (shown in Fig. |^). Whilst the 
magnetic structure of Mnl and Mn2 is fixed, the spin di¬ 
rections of Mn3 ions are allowed to fully relax. The above 
scheme is employed for all calculations unless otherwise 
stated. 

It is well known that both magnetic and electric prop¬ 
erties from GGA-\-U calculations are very sensitive to 
the Hubbard parameter U. Experiments show that 
CaMn 70 i 2 is a bad insulator at room temperature with a 
band gap <0.38 eV. In our calculations, we fixed J = 1 
eV and tried several different U values. We find that 
for <3 eV, the experimental spin configuration be¬ 
comes metallic, while U = 3 eV produces a small band 
gap ~ 0.15 eV. There is no other qualitative difference 
among the results calculated with different U values, so 
here, unless otherwise stated, we show the results calcu¬ 
lated with U = 3 eV and J = 1 eV. 


III. STRUCTURAL MODULATION AND 
ORBITAL HELIX 

At 150K, the modulated structure can be character¬ 
ized by a centrosymmetric four-dimensional space group 
R 3 ( 007 ) with propagation vector Qc = (0, 0, 2.077).I^This 
structural modulation is dominated by variations of Mn2- 
O bonds, with smaller variations of Mnl-0 and Mn3-0 
bonds. The MnOe octahedral crystal field splits Mn 3d 
energy levels into higher Cg and lower t 2 g states. Since 
Mn2 has a 3d^ configuration (Mn^+), the t 2 g states are 
fully occupied while the eg states are partially occupied 
with parallel spins. The structural modulation leads 
to local deformations of the oxygen octahedra, which 
can be characterized by two distortion modes: tetrag¬ 
onal (Qs) and orthorhombic (Q 2 )- Qs and Q 2 can 
be further expressed as Qs = {2l\/Q)(2Z — X — Y), 
Q 2 = {2/y/2){X — T), where X, Y, and Z correspond 
to Mn-0 bond lengths in the x,^and 2 : directions de¬ 
fined in the local fram^ (see Fig. |h. The general distor¬ 
tion can then be written as jQ) = cos^lQs) + sin a | ( 32 ), 
where a is a mixing angle. According to the Jahn- 
Teller effect, the wave function of an Cg state is given 
by |4>) = cos{a/2)\3z^ — r^) -h sin((a/2)|x^ — where 

tan((a) = Q 2 lQ?> = (^^-x-y) ^ show that 

a = 27t/3 and 47 r /3 correspond to the occupation of 
3^2 _ ^2 ^y 2 _ ^2 orbitals respectively. 

In the average (unmodulated) structure, the Mn^+Oe 
octahedra of the Mn2 site are compressed, giving two 
short bonds along 2 : and four long bonds in the xy plane. 
At variance with this situation, the modulated structure 
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FIG. 1: Experimental bond-length modulation in CMO : 
Mn2-0 bond lengths along the local octahedral x, y and 2 ; 
(see inset) as a function of plotted across 15 unit cells 
along the c axis. 



FIG. 2: Partial density of states for the Mn2 sites of structures 
Al, A2, and A3. The Fermi energy is at Ef = 0. 


leads to an incommensurate antiphase variation in the 
bond lengths in the xy plane, consistent with a rotation 
of a and therefore a periodic modulation of the occupa¬ 
tion of the 3x‘^ — and 3y‘^—r‘^ orbitals. The modulation 
of the Mn2-0 bond lengths is shown in Fig. Since the 
propagation vector is incommensurate, it is not possi¬ 
ble to study the full structural modulation directly from 
first-principles calculations, since the corresponding su¬ 
per cell is infinitely long and any realistic commensurate 
approximation would be computationally prohibitive. As 
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FIG. 3: (a) Side view of one unit cell of average crystal struc¬ 
ture and magnetic configuration. Yellow (large) and blue 
(small) spheres represent Ca and Mn3 ions respectively. Red 
(dark) and green (light) spheres with arrows represent Mnl 
and Mn 2 ions respectively. Arrows denote magnetic moments 
of Mn ions. Nearest neighbour exchange interations are indi¬ 
cated by dotted lines, (b) Top view and (c) side view of one 
cluster centered at Mn32 atoms. Dotted lines indicate the 
splitting of J 3 and J 4 in the modulated structure. 


an alternative, we study a series of models designed to 
capture the local electronic structure in the presence of 
the modulation. In this approach, an individual unit cell 
along the c direction is selected and then built up as a 
crystal using periodic boundary conditions. The experi¬ 
mental structural parameters are used, unless otherwise 
stated. Orbital order is then investigated by determin¬ 
ing the orbital occupation of the Mn2 atom centred at 
( 1 / 2 , 1 / 2 , 1 / 2 ), so as to minimise the influence of the ar¬ 
tificial unit cell boundary. In order to keep the com¬ 
putational load tractable, we choose to perform calcula¬ 
tions on three typical unit cells marked with Rz = 0,3,6 
in Fig. denoted here as Al, A 2 and A3, respectively. 
These points correspond to the two extrema of the mod¬ 
ulation (Al and A3) and the intermediate nodal point 
(A 2 ). The a values of Ai, A 2 , and A 3 structures are 
128°, 183° and 223°, which correspond to Cg occupation 
with dominant 3x^ — r^, and 3y^ — characters 

respectively. A standard way to understand the local or¬ 
bital occupation is through calculation of the partial den¬ 
sity of states (PDOS). The calculated PDOS are shown 
in Fig.[^ From Fig.[^ we can see that the Cg orbitals are 
separated at the Fermi level by a small gap. The respec¬ 
tive Cg orbitals of Al, A 2 , and A3 clearly have a domi¬ 
nant 3x‘^ — r^, x‘^ —y‘^ and 3y‘^ — occupation below the 
energy gap, respectively, which is in excellent agreement 
with predictions based on the Jahn-Teller mechanism. 


IV. MAGNETIC STRUCTURE 

In addition to the pure electronic structure, magnetic 
ordering in CMO is also expected to be strongly affected 
by the structural modulation. To understand this as¬ 
pect, we performed additional calculations on the series 
of structure models described above, investigating specif¬ 
ically the coupling between structure and magnetism. In 
each unit cell of CaMn 70 i 2 , there are a total of three Mn3 
atoms located in different layers, which we denote here 
as Mn31, Mn32, and Mn33, respectively, as shown in Fig. 
I^a). Each Mn3 ion forms a cluster with 6 nearest neigh¬ 
bour (NN) Mnl and 6 NN Mn2, as illustrated in Fig. 
I^b) and (c). The directions of the magnetic moments 
of Mn31, Mn32, Mn33 with respect to the neighbouring 
Mnl and Mn2 layers are expressed through the angles 
6 >i, 6 > 2 , and 6 ^ 3 . 

We employed a Heisenberg-like model including 
both NN symmetric exchange (SE) interactions and 
antisymmet ric Dzy aloshinskii-Moriya (DM) exchange 
inter act ions HI] xhe magnetic interaction associated 

with one pair of Mn ions joined by a common oxygen 
ligand ion can be written as 

H = JijSi ■ Sj + 'Dij ■ (Si xSj), (1) 

where the first term is the Heisenberg exchange energy 
and the second term is the antisymmetric DM exchange 
energy which originates from spin orbit coupling (SOC). 
Si and Sj are spins of Mn^, Muj, D^j = /?Ui X Uj is the 
DM vector, where and Uj are vectors connecting the 
two Mn sites with the oxygen ligand, and ^ is a. cou¬ 
pling constant. Eor NN SE interactions, we use the same 
notations as in Ref. m As shown in Eig. there 
are 7 inequivalent exchange interactions in the average 
structure Ji — Jj, among which J 3 and J 4 correspond 
to Mn3-Mnl and Mn3-Mn2 super-exchange interactions, 
respectively. In the modulated crystal structure, the ex¬ 
change interactions are modulated along the c direction. 
The resulting splitting of J 3 and J 4 is illustrated for the 
Mn32 cluster in Eig. |^c) where in one cluster J 3 and 
J 4 are relabelled J 31 , J 41 , and J 32 , J 42 , corresponding to 
interactions with the upper and lower Mnl and Mn2 lay¬ 
ers, respectively. The DM vectors D are denoted using 
the same subscripts as the SE interactions. According to 
Eq. 0 ’ the exchange energy of the Mn32 cluster with 
the empirical spin directions fixed for Mnl and Mn2 can 
be written as 


E = 


Ja COS (I - ^2) + Jb COS (^ + 6*2) 
+Da sin (I - 6 * 2 ) - Db sin (| + 6 * 2 ) 

cos(02) 


l{Ja + Jb) + ^{Da-Db) 


+ 


^{Ja-Jb)-l{Da+Db) 


sin( 6 » 2 ). ( 2 ) 
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where Ja = 3(J31+J42), Jh = 3(J32 + J4l): 

1 ) 42 ), and Di) = 3(1^32 +4^|i). By minimizing the energy 
with respect to O 2 , we obtain 


tan( 6 > 2 ) 


V3{Ja - Jh) - {Dg + 4 ^ 5 ) 

+ Jh) + V^{Da — Di)) 


(3) 


These equations can be generalized to other Mn3 clusters 
in the whole lattice. In the average structure 0i = O 2 = 
O 3 = 0 is enforced by symmetry, and Mn3 atoms occupy 
sites with 3 point symmetry, which implies that Ja = Jh 
and ]^ = Di). In the absence of SOC, Da = Dij = ^ and 
Eq. ( 3 I becomes tan(^ 2 ) = ; from which we can 

easily obtain 6 > = 0. On the other hand, including SOC 
will cause a deviation from 0° due to DM interactions. 
Our ab initio calculations of the average structure give 
0 = 0 without SOC and 0 = —4.1° when SOC is taken 
into account. This supports the validity of the above 
analysis. 

As already stated, in the fully modulated atomic struc¬ 
ture, the Mn3-Mnl and Mn3-Mn2 bonds are modu¬ 
lated along the c direction. As in Sec. |ml we study 
the modulation of the Mn3 spin direction by perform¬ 
ing bulk calculations on periodic structures constructed 
from unit cells extracted from the modulated structure. 
Once again, to reduce the impact from artificial neigh¬ 
bouring cells along c direction, we only investigate the 
Mn32 atom located in the middle layer of each unit cell 
with 2 ; = 1/2. To illustrate the structural modulation 
around the Mn32 site quantitatively, Mn32-Mn2 bond 
lengths and Mn32-0-Mn2 bond angles are plotted in 
Fig. |^a)(b). It can be seen that both bond lengths 
and bond angles vary with the wave vector qc, as ex¬ 
pected. The modulation of Mn32-Mnl has a very similar 
behavior. We can therefore naturally assume that both 
Ja and J 5 — see Eq. Q — will be modulated with the 
wave vector qc, i.e. JcJv) = Jq + AJsin(qc • r + 0 + 7 /), 
J 5 (r) = Jo+AJ sin(qc'r-F 0 ), where AJ is the magnitude 
of the change in J as a result of the modulated structure, 
0 is a phase shift and 7 / is the phase difference between 
these two modulations. We find that, in the modulated 
structure, the vectors D/^d connecting the Mn3 layers to 
the NN Mn2 and Mnl layers are only slightly modulated 
and always have the same sign, indicating a weak mod¬ 
ulation of both Da and D 5 for all three Mn3 sites. The 
ab-initio calculated O 2 are shown in Fig. Qc). We can 
see that the O 2 modulation is essentially identical in the 
presence and in the absence of SOC, indicating that the 
dominant factor at play is the effect of bond lengths and 
angles on the Heisenberg exchange. To compare with 
the average structure, we can consider for example the 
A 3 structure {z = 6 ). For A 3 , Ja = Jb, as in the aver¬ 
age structure, and hence ^2 ~ 0 without SOC, as before. 
The only effect of the SOC is to introduce a constant 
O 2 offset of approximately —4°, consistent with the opti¬ 
mized 0 = —4.1° in the average structure. In order for 
the offset to be constant along c, Eq. (§ imposes that 
Da and D 5 should also be weakly modulated along the 


structural modulation. Consequently, we can safely set 
Da = D}) = D and re-write Eq. (§ as 


tan ^ 2 ( 1 *) 


V3[Ja{r)-Mr)]-2D 

Ja(r) + J(,(r) 


(4) 


By fitting Eq. 0 to the ab initio values of 02 in Fig.|^c), 
we obtain AJ/Jq = 0.49, 0 = 0.407r, 7 / = 0.357r, 
AD/Jq = 0.08. The quality of this fit is excellent as 
can be seen in Fig. Qc). The modulation of J is about 
half of Jo, preserving the sign of the exchange interac¬ 
tions throughout the whole crystal. To double check the 
validity of the above analysis, we calculated the mag¬ 
netic moment magnitudes of Mn32 along the modula¬ 
tion, which are shown in Fig. Qd). It can be seen that 
the magnetic moments are around 2.18 i^b with a slight 
modulation of amplitude ^ 0 . 1 //^ 5 which agrees very well 
with experiments.^ 

One important consequence of the O 2 modulation is to 
introduce an additional magnetic propagation vector for 
the Mn3 sublattice. A Fourier analysis of the Mn3 spin 
structure obtained from Eq. 0 gives three modulation 
vectors q^i, qc + qmi and qc - qmi- Since qc = 2 q^i, 
Qc “ Qm! = Qm! SO that the first and third propagation 
vectors are equal, whereas qc + qmi = 3qmi- This leads 
us to predict that the magnetic structure of the Mn3 
sublattice in the AFMl phase is anharmonically modu¬ 
lated with two propagation vectors, qmi and 3qmi,due 
to orbitally-mediated coupling with the structural mod¬ 
ulation. In this context, the term “anharmonic” refers 
purely to the presence of higher-order harmonics of the 
magnetic propagation vector qmi- As we have seen, this 
conclusion is independent on whether or not the SOC is 
taken into account, since this does not change the peri¬ 
odicity of the O 2 modulation. 


V. ELECTRIC POLARIZATION 


In this section, we describe our calculations of the elec¬ 
tric polarization in the AFMl phase of CaMn 70 i 2 per¬ 


formed with the same parameter setup as in Section IV 


We initially performed calculations on the average atomic 
structure for comparison, where the magnetic structure 
is only characterized by the spin direction 0 of Mn3, as 
before. The electric polarization calculated as a function 
of the angle 0 is shown in Fig. Qa). The distinct contri¬ 
butions to the polarization from the symmetric exchange 
striction (Pes) and inverse DM effects (Pdm) are iden¬ 
tified by turning on and off the SOC. Fig. Qa) shows 
that Pes — 6 >, while Pdm is approximately constant. 
As already remarked,!^ Pes must be zero for 0 = 0 be¬ 
cause magnetic structures with opposite chiralities are 
related by a global rotation in spin space. These calcu¬ 
lated results can be explained using the analytical model 
described in Section |IV| as follows. By considering the 
same Mn32 cluster illustrated in Fig. Qc), the contribu¬ 
tion from symmetric exchange striction can be written 
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FIG. 4: (a)(b) Experimental data of modulated Mn32-Mn2 
bond lengths and bond angles along the c direction.^ (c) O 2 
from ah initio calculations. Red and black lines are htted to 
Eq. Dashed line is drawn at ^2 = 0. (d) The correspond¬ 
ing magnetic moments of the Mn3 ions. 
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EIG. 5: Electric polarization of (a) the average structure and 
(b) the GM structure. Symbols denote results obtained from 
ab initio calculations. Lines represent hts to Eq. § and 

Eq. (|^. 


spectively. This contribution can therefore be written as 

ftM(^2) / 341 S 0 X (Si + S 2 + S 3 ) 

+ ^42^0 X (S 4 + S 5 + Se) 

+ ^3lSo X (S7 + S8 + S9) 

+ ^ 32^0 X (Sio + Sii + S12) 

= Pa sin - $ 2 ^ + Pb sin 


TT + ^2 


= +/36)cOs(6»2) 

+ ^(/3a - ^ 6 )sin( 6 » 2 ), 


(6) 


-Pes(^* 2 ) 


oc 741S0 ■ (Si + S2 + S3) 

- 742S0 ■ (S4 + Ss + Se) 

+ 731S0 • (S7 + 83 + Sg) 

~ 732S0 • (Sio + Sii + S12) 

= 7 aCos(| - 6»2) + 7 bCOs(| +6*2) 

= ^(7a +76 )cOs(6»2) 

a/3 

+ ^(7a-76)sin(02), 


( 5 ) 


where 731,732,7415 and 742 are exchange striction pa¬ 
rameters. In the linear approximation, the parameters 
are proportional to the first-order derivative of the cor¬ 
responding J’s with respect to the ionic displacements.^ 
In Eq. we have defined 7^ = 3(731 — 742), and 
76 = 3(741 — 732). Likewise, the contribution from SOC 
can be accurately described as a DM striction effect, i.e. 
due to the imbalance in the DM interactions between 
Mn 3 and the upper and lower Mnl and Mn 2 layers, re- 


where Psi^ P32^ Pai^ and P42 are coupling constants corre¬ 
sponding to D31, D32, D41, and D42, respectively, and we 
have defined pa = 8(^31 + ^42), and Pb = 8(^41 + ^32)- 
The contribution to the total polarization from Mn 31 
and Mn 33 clusters could be derived in the same way. 
The total polarization is then Ptot = ^es + Pbm , where 
Pes = PEsi^i) and Pdm = ft>M(^i)- 

In the average structure 0 i = O2 = O3 = 0 , ^a = ~lb 
and Pa = Pbi hence Eq. (§ may be simplified to 
Pes C)C sin(6>), and Eq. ^ to Pdm cos{ 0). The corre¬ 
sponding fits are also drawn in Fig. [^as lines for com¬ 
parison. We can see that there is excellent agreement 
between the analytical expression and the results calcu¬ 
lated ab initio. However, it is worth noting that the total 
polarization calculated for the optimized 0 = 4 . 1 ° is only 
- 0.46 mC/m^, about five times smaller than the experi¬ 
mental value. We will come back to this after examining 
the electric polarization in a modulated structure. 

As the electric polarization is a bulk property, we can¬ 
not approximate the full modulated structure by repli¬ 
cating a single unit cells of interest, as we have done 
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FIG. 6: Schematic illustrations of the splitting of the J 4 
exchange interactions in the CM structure. The color scheme 
is the same as in Fig. 


in the previous sections. We therefore investigate a 
commensurately modulated (CM) atomic structure with 
qj, = ( 0 , 0 , 2 ), constructed from the experimental data.^^ 
It is worth emphasising the subtly different approach we 
take here: in sections III and IV, we “slice off” a sec¬ 
tion of the incommensuratey modulated structure corre¬ 
sponding to a hexagonal conventional unit cell, which is 
a supercell of the rhombohedral unit cell, and replicate 
it. This results in locally exact bond lengths and angles, 
but also in a slight mismatch at the boundary between 
unit cells. In this section, we fix the 4-dimensional frac¬ 
tional coordinates of the modulation to the experimental 
values, but we approximate the propagation vector with 
the commensurate value qj, = (0,0,2). This also results 
in a hexagonal supercell of the average structure. The lo¬ 
cal atomic positions, and consequently bond lengths and 
angles, are everywhere only an approximation of the real 
structure, but the modulation is smooth and there is no 
discontinuity at the unit cell boundary. The CM atomic 
structure is constructed based on the cell at Rz = 0 in 
Fig.[T} In the super-cell of this CM structure, each type 
of exchange interaction is split into three distinct interac¬ 
tions, i.e. Ji splits into J^i, 7 ^ 2 , Jis- By way of example, 
the splitting of J 4 is illustrated in Fig. Inversion sym¬ 
metry in the CM structure is still preserved with Mn32 
occupying the inversion centre, while Mn31 and Mn33 
occupy two sites related by inversion symmetry. 

According to the symmetry of the positions occupied 
by Mn3 atoms, we have 6^2 = 0 and 0i = —63 for mag¬ 
netic structures without SOC. Indeed, our calculations 
give values of 0i = 25°, O 2 = 0, 6^3 = —25°. By in¬ 
cluding SOC we find 0i = 20 °, O 2 = —5°, O 3 = —30°, 
where all the three 0 values show a constant rotational 
offset of = —5° compared to that calculated without 
SOC, consistent with our analysis and our calculated re¬ 
sults in Sec. |lv| In the CM structure, a straightforward 
calculation according to Eq. <§ and Eq. shows that 
Pes oc sin(A^) and Pdm oc cos(A^). To further study 
this relation, the polarization calculated by fixing vari¬ 


ous A 6 > are shown in Eig. [^(b). Here, a perfect fit to the 
analytical sine and cosine function can be obtained. The 
contributions from Mn31 and Mn33 to Pes partly cancel 
each other due to their opposite rotation directions. The 
overall behaviour of the electric polarization is then very 
similar to that of the average structure, except for the 
fact that 0 of the averaged structure is replaced by the 
global rotational offset A^. In both cases, the total polar¬ 
isation results from fine tuning of Pes and Pdm • For the 
ab initio optimized AO = —5°, we calculate Pes = +2.03 
mC/m^ and Pdm = —2.07 mC/m^, resulting an almost 
exact cancellation of the total polarization (Ptot = —0.04 
mC/m^ — two orders of magnitude smaller than the ex¬ 
perimental value). We also performed calculations with 
CM atomic structures constructed from other cells of the 
ICM structure. Polarizations from different CM struc¬ 
tures show similar modulation with Rz as in Eig.[^ while 
the total polarization calculated by averaging the polar¬ 
izations from each CM structure is consistent with that 
obtained using the CM cell from Rz 

One probable reason for this coincidental cancellation 
is that in our calculation the atomic structure is only ap¬ 
proximated, as explained above, and in not fully relaxed. 
We therefore repeated our ab-intio calculations with re¬ 
laxed atomic positions of the CM structure but keeping 
the experimental lattice parameters. After relaxation, 
AO slightly increases to —5.6° and the total polarization 
is greatly enhanced to Ptot = —1.19 mC/m^, much closer 
to the experimental value. It is worth reminding that the 
sign agreement between measured and calculated polar¬ 
izations is extremely difficult to determine, since it re¬ 
quires establishing experimentally the absolute chirality 
of the magnetic structure. The atomic relaxation changes 
the Mn-Mn bonds and hence the corresponding coupling 
parameters 7 and leading to an increase in the con¬ 
tribution from the DM interaction relative to that from 
exchange striction. Another possibility is that the cal¬ 
culated value of AO is too small, possibly due to an 
underestimation of the DM interaction. If we enforce 
AO = — 10 °, the contribution from exchange striction 
prevails and a total polarization of Ptot = 2.4 mC/m^ is 
obtained, which is in good agreement with experiments^. 
Erom this analysis, it is clear that the DM interaction 
plays a key role in producing macroscopic electric polar¬ 
ization. In the absence of DM interaction, AO = 0, and 
both Pes and Ptot vanish. When the DM interaction is 
present, a net electric polarization arises with contribu¬ 
tions from both Pes and Pdm- This conclusion, based on 
the analysis of our calculated data, is robust and largely 
insensitive to the calculation details. 


VI. CONCLUSIONS 

We have presented a series of ab initio calculations 
on multiferroic CaMn 70 i 2 , based on the experimentally 
determined modulated crystal structure and designed to 
provide insight into the electronic and magnetic struc- 
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tures and the nature of the magnetically induced elec¬ 
trical polarization of this material. We also presented a 
series of analytical models expressed in terms of a small 
number of parameters, which provide excellent fits to the 
ab-intio calculations. By investigating the 3d orbital oc¬ 
cupation using a PDOS analysis, we confirmed the 3d 
Cg orbital modulation of Mn2 atom previously proposed. 
We further investigated the detailed magnetic structure 
of the Mn3 (Mn^+) sublattice in the multiferroic AFMl 
phase and its interactions with the Mn^+-containing Mnl 
and Mn2 sublattices. We find that the Mn3 magnetic 
structure is anharmonically modulated due to orbitally- 
mediated coupling with the structural modulation — a 
key ingredient in stabilising the observed magneto-orbital 
helices. Finally, we studied the magnetoelectric polariza¬ 
tion both in the average (unmodulated) structure and in 


a commensurate approximation of the modulated struc¬ 
ture. The Dzyaloshinskii-Moriya interaction was found 
to play a key role in producing magnetoelectric polariza¬ 
tion, while the total polarisation results from a delicate 
interplay of inverse DM and exchange-striction contribu¬ 
tions. 
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